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We study the gauge/gravity dual of a finite temperature field theory at finite isospin chem- 
ical potential by considering a probe of two coincident D7-branes embedded in the AdS- 
Schwarzschild black hole background. The isospin chemical potential is obtained by giving 
a vev to the time component of the non-Abelian gauge field on the brane. The fluctuations 
of the non-Abelian gauge field on the brane are dual to the SU{2) flavor current in the field 
theory. For the embedding corresponding to vanishing quark mass, we calculate all Green 
functions corresponding to the components of the flavor current correlator. We discuss the 
physical properties of these Green functions, which go beyond linear response theory. In par- 
ticular, we show that the isospin chemical potential leads to a frequency-dependent isospin 
diffusion coefficient. 
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I. INTRODUCTION 

Over the past years, there have been a number of lines of investigation for describing QCD- 
like theories with gravity duals. In this way, considerable progress towards a gauge/gravity dual 
description of phenomenologically relevant models has been made. One of these lines of investi- 
gation is the gravity dual description of the quark-gluon plasma obtained by applying AdS/CFT 
to relativistic hydrodynamics |QI Q, 0, 3, H, Is 0]. The central result of this approach is the cal- 
culation of the shear viscosity from AdS/CFT. More recently, an R charge chemicalpotential has 
been introduced by considering gravity backgrounds with R charged black holes [@,& 10], and 



also the heat conductivity has been calculated by considering the R current correlators in these 
backgrounds |8|]. 

A further approach to generalizing the AdS/CFT correspondence to more realistic field theories 



is the addition of flavor to gravity duals via the addition of probe branes [|llLll2lll3Lll4l ll5ll. This 
allows in particular for the calculation of meson masses. 

These two approaches have been combined in order to study the flavor contribution to finite 
temperature field theories from the gravity dual perspective. This began with lfl3ll where the 
embedding of a D7-brane probe into the AdS-Schwarzschild black hole background was studied 
and a novel phase transition was found, which occurs when the D7 probe reaches the black hole 
horizon. This transition was shown to be of first order in lllql (see 11141 for a similar transition 



in the D4/D6 system), and studied in further detail in 11171 11811 . Related phase transitions appear 



in 



191. 201. 1211. Mesons in gravity duals of finite temperature field theories have been studied in 

BE3MI3. 



Recently, in view of adding flavor to the quark-gluon plasma, the flavor contribution to the 



shear viscosity has been calculated in Il26l l27h. where it was found that r/f un d oc \N c NfT 3 . 

For a thermodynamical approach in the grand canonical ensemble, the inclusion of a chemical 
potential and a finite number density is essential. In rf28ll . an isospin chemical potential was in- 
troduced by considering two coincident D7 probes, and by giving a vev to the time component of 
the SU (2) gauge field on this probe. This was shown to give rise to a thermodynamical instability 
comparable to Bose-Einstein condensation, compatible with the field-theoretical results of ^ 
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For the gauge/gravity dual analysis, a potential generated by an SU{2) instanton on the D7 probe 
in the gravity background was used 



A baryon chemical potential (j,r is obtained by turning on the diagonal U(l) C U(Nf) gauge 
field on the D7-brane probe 113411 . Contributions to the D7-brane action arise from the deriva- 
tive of this U(l) gauge field with respect to the radial direction. The effects of this potential on 



the first-order phase transition described above have been studied in ||34fl . where regions of ther- 
modynamical instability have been found in the (T, phase diagram. - For the DA/ D8/ D8 
Sakai-Sugimoto model [15, [lH], the phase transitions in presence of a baryon number chemical 



potential, as well as physical processes such as photoemission and vector meson screening, have 
been studied in BSij]. 



A related approach has been used to calculate the rate of energy loss of a heavy quark moving 
through a supersymmetric Yang-Mills plasma at large coupling [39]. In this approach the heavy 
quark is given by a classical string attached to the D7-brane probe. - A first study of flavors in 
thermal AdS/CFT beyond the quenched approximation, i.e. with Nf ~ N c , was performed in ll40ll . 

Here we study finite-temperature field theories with finite isospin chemical potential by consid- 
ering two coincident D7-brane probes in the Lorentzian signature AdS-Schwarzschild black hole 



background. As in II28I1 . we introduce an isospin chemical potential by defining 



for the time component of the SU(2) gauge field on the two coincident D7-branes. This constant 
chemical potential is a solution to the D7-brane equations of motion and is present even for the 
D7-brane embedding corresponding to massless quarks. We consider small fi, such that the Bose- 
Einstein instability mentioned above, which is of order C(/i 2 ), does not affect our discussion here. 

For simplicity we consider only the D7 probe embedding for which the quark mass vanishes, 
m = 0. This embedding is constant and terminates at the horizon. We establish the SU(2) non- 
Abelian action for a probe of two coincident D7-branes. We obtain the equations of motion for 
fluctuations about the background given by (11.11) . These are dual to the SU{2) flavor current J^ a . 
We find an ansatz for decoupling the equations of motion for the different Lorentz and flavor com- 
ponents, and solve them by adapting the method developed in [0, [3[] . This involves Fourier trans- 
forming to momentum space, and using a power expansion ansatz for the equations of motion. We 
discuss the approximation necessary for an analytical solution, which amounts to considering fre- 
quencies with cu < fi < T. With this approach we obtain the complete current-current correlator. 
The key point is that the constant chemical potential effectively replaces a time derivative in the ac- 
tion and in the equations of motion. In the Fourier transformed picture, this leads to a square-root 
dependence of physical observables on the frequency, y/uj. This non-linear behavior goes beyond 
linear response theory. We discuss the physical properties of the Green functions contributing to 
the current-current correlator. In particular, for small frequencies we find a frequency-dependent 
diffusion coefficient D(co) oc ^^fojjji. Whereas frequency-dependent diffusion has - to our 
knowledge - not yet been discussed in the context of the quark-gluon plasma, it is well-known in 
the theory of quantum liquids. For instance, for small frequencies the square-root behavior we find 
agrees qualitatively with the results of ll4lL l42ll for liquid para-hydrogen. Generally, frequency- 
dependent diffusion leads to a non-exponential decay of time-dependent fluctuations, as discussed 



for a classical fluid in 114311 



Physically, the isospin chemical potential corresponds to the energy necessary to inverting the 
isospin of a given particle. Within nuclear physics, such a chemical potential is of relevance 
for the description neutron stars. Moreover, isospin diffusion has been measured in heavy ion 
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reactions 11441. 14511 . - For two-flavor QCD, effects of a finite isospin chemical potential have been 
discussed for instance in [|46l 1471. 14811 . The phase diagrams discussed there are beyond the scope 
of the present paper. We expect to return to similar diagrams in the gauge/gravity dual context in 
the future. 

This paper is organized as follows. In section 2 we summarize the AdS/CFT hydrodynamics 
approach to calculating Green functions, which we use in the subsequent. Moreover we comment 
on frequency-dependent diffusion within hydrodynamics. In section 3 we establish the D7 probe 
action in presence of the isospin chemical potential, derive the corresponding equations of motion 
and solve them. In section 4 we obtain the associated Green functions in the hydrodynamical 
approximation. We discuss their pole structure and obtain the frequency-dependent diffusion co- 
efficient. We conclude in section 5 with an interpretation of our results. An explanation of our 
notation as well as a series of calculations are relegated to a number of appendices. 



H. HYDRODYNAMICS AND ADS/CFT 

Thermal Green functions have proven to be a useful tool for analyzing the structure of hydro- 
dynamic theories and for calculating hydrodynamic quantities such as transport coefficients. For 
instance, given a retarded current correlation function G(k)^ u in Minkowski space, the spectral 
function can be written in terms of its imaginary part, 

X^ik) = -2lmG^(k). (2.1) 



For the gravity dual approach, this is discussed for instance in lf7L|49H. 

In this paper we use the gauge/gravity dual prescription of |3|] for calculating Green functions 
in Minkowski spacetime. For further reference, we outline this prescription in the subsequent. It is 
based on the AdS/CFT-correspondence relating supergravity fields A in a black hole background 
to operators J in the dual gauge theory. The black hole background is asymptotically Anti-de 
Sitter space and places the dual field theory at finite temperature. This temperature corresponds 
to the Hawking temperature of the black hole, or more generally speaking, of the black branes. 
Starting out from a classical supergravity action S c \ for the gauge field A, according to IH] we 
extract the function B{u) (containing metric factors and the metric determinant) in front of the 
kinetic term (d u A) 2 , 

S cl = J dud 4 xB(u) {d u Af + ... (2.2) 

Then we perform a Fourier transformation and solve the linearized equations of motion in momen- 
tum space. This is a second order differential equation, so we have to fix two boundary conditions. 
The first one at the boundary of AdS at u = can be written as 

A(u,k) = f(u,k)A M y(k), (2.3) 

where A hAy {k) is the value of the supergravity field at the boundary of AdS depending only on 
the four flat boundary coordinates. Thus by definition we have f(u,k) =1. For the other 
boundary, located at the horizon u = 1, we impose the incoming wave condition. This requires 
that any Fourier mode A(k) with timelike k can travel into the black hole, but is not allowed to 
cross the horizon in the opposite direction. For spacelike k, the components of A have to be regular 
at the horizon. Then the retarded thermal Green function is given by 



G(u,q) = -2B(u)f(u,-k)d u f(u,k) 



(2.4) 

u-s-0 
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The thermal correlators obtained in this way display hydrodynamic properties, such as poles 
located at complex frequencies. Generically, for the R current component correlation functions 
calculated from supergravity, there are retarded contributions of the form 

G(u,q)oc- 1 — 2 . (2.5) 

ioj — Dq 2 

This may be identified with the the Green function for the hydrodynamic diffusion equation 

d J (t,x) = DV 2 J (t,x), (2.6) 

with J the time component of a diffusive current. D is the diffusion constant. In Fourier space 
this equation reads 

iuJ (uj, k) = Dk 2 J (u, k) . (2.7) 

In position space, this corresponds to an exponential decay of J with time. 

For the non-Abelian case with an isospin chemical potential, in sections [III] and [IV] we will 
obtain retarded Green functions of the form 

G{u>,q)<x- ^-j. (2.8) 

Retarded Green functions of this type have been discussed for instance in ll43ll . (|2.8I) corresponds 
to frequency-dependent diffusion with coefficient D{uj), such that (|2.7I) becomes 

iujJ {u, k) = D{u)k 2 J (u, k) . (2.9) 

In our case, Jo is the averaged isospin at a given point in the liquid. 

This is a non-linear behavior which goes beyond linear response theory. In particular, when 
Fourier-transforming back to position space, we have to use the convolution for the product D ■ J 
and obtain 

t 

d J (t,x) + V 2 J dsJ {s,x)D(t- s) = (2.10) 

— oo 

for the redarded Green function. This implies together with the continuity equation d J + V- J = 
0, with J the three-vector current associated to J , that 

J = _V( J D*J ), (2.11) 

where * denotes the convolution. This replaces the linear response theory constitutive equation 
J = -DVJq. Note that for D(t - s) = D5(t - s) with D constant, ([2.10b reduces again to (f2T6l) . 



III. SUPERGRAVITY BACKGROUND AND ACTION 

A. Finite temperature background and brane configuration 

We consider an asymptotically AdS 5 x S 5 spacetime as the near horizon limit of a stack of N c 
coincident D3-branes. More precisely, as in [2D, our background is an AdS black hole, which is 
the geometry dual to a field theory at finite temperature. The Minkowski signature background is 

ds 2 = dx 2 + dx 2 + dx 2 dx 2J du 2 R 2 dQ 2 

u v 3y Au 2 f(u) ~ 

b*R* 1 } 



< u < 1, Xi e R, C, 



0123 



u 



2 
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with the metric dfi 2 , of the unit 5-sphere, and the function f(u), AdS radius R and temperature 
parameter b given in terms of the string coupling g s , temperature T, inverse string tension a' and 
number of colors iV c by 

f(u) = 1 - u 2 , R A = 4irg s N c a' 2 , b = irT. (3.2) 

The geometry is asymptotically Ad S 5 x S 5 with the boundary of the AdS part located at u = 0. 
At the black hole horizon the radial coordinate u has the value u — 1. 

Into this ten-dimensional spacetime we embed Nf = 2 coinciding D7-branes, hosting flavor 
gauge fields A^. The embedding we choose extends the D7-branes in all directions of AdS space 
and wraps an S 3 on the S 5 . In this work we restrict ourselves to the most straightforward case, 
that is the embedding of the branes through the origin along the AdS radial coordinate u. This 
corresponds to massless quarks in the dual field theory. On the brane, the metric in this case simply 
reduces to 

b 2 R 2 f? 2 

d, 2 = (-/(«) dxl + dxl + dxj + dxl) + du 2 + RW 3 , 

u Au z j{u) (3.3) 

< u < 1, Xi E R. 

Due the choice of our gauge field in the next subsection, the remaining three-sphere in this metric 
will not play a prominent role. 

The table below gives an overview of the indices we use to refer to certain directions and 
subspaces. 





AdS 5 S 3 


coord, names 


n 1 2 ^ 
X X x x u - 




/i, v.. . 


indices 





6. Introducing a non-Abelian chemical potential 

A gravity dual description of a chemical potential amounts to a non-dynamical time component 
of the gauge field A^ in the action for the D7-brane probe embedded into the background given 
above. There are essentially two different ways to realize a non-vanishing contribution from a 
chemical potential to the field strength tensor F = 2d[ fM A u ] + f abc A b ^Al. The first is to consider 
a M-dependent baryon chemical potential for a single brane probe. The second, which we pursue 
here, is to consider a constant isospin chemical potential. This requires a non-Abelian probe brane 
action and thus a probe of at least two coincident D7-branes, as suggested in Il28ll . Here the time 
component of the gauge field is taken to be 

A = A a t a , (3.4) 

where we sum over indices which occur twice in a term and denote the gauge group generators by 
t a . The brane configuration described above corresponds to an SU(Nf) gauge group with Nf = 2 
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on the brane, which corresponds to a global SU (Nf) in the dual field theory. For Nf = 2, the 
generators of the gauge group on the brane are given by t a = with Pauli matrices a a . We will 
see that (|3.4I) indeed produces non-trivial new contributions to the action. 

Using the standard background field method of quantum field theory, we will consider the 
chemical potential as a fixed background and let the gauge fields fluctuate. We single out a par- 
ticular direction in flavor space by taking Aq = fi as the only non-vanishing component of the 
background field. From now on we use the symbol A a v to refer to gauge field fluctuations around 
the fixed background, 

A a u ->fi5 u0 5 a3 + A a u . (3.5) 

We pick the gauge in which A u = and assume that A^ = for /i = 5, 6, 7. Due to the 
symmetries of the background, we effectively examine gauge field fluctuations A^ living in the 
five-dimensional subspace on the brane spanned by the coordinates xq, x±, x-i, £3 and by the radial 
AdS coordinate u. The magnitude of all components of A and the background chemical potential 
/i are considered to be small. This allows us to simplify certain expressions by dropping terms of 
higher order in A and in the chemical potential ji. 



C. Dirac-Born-Infeld action 

The action describing the dynamics of the flavor gauge fields in the setup of this work is the 
Dirac-Born-Infeld action. There are no contributions from the Chern-Simons action, which would 
require non-zero gauge field components in all of the 4,5,6,7-directions. As mentioned, we con- 
sider the D7 probe embedding whose asymptotic value at the boundary is chosen such that it cor- 
responds to vanishing quark mass, m = 0. The metric on the brane is then given by (13.31) . Since 
we are interested in two-point correlators only, it is sufficient to consider the action to second order 



in a', 



Sdi = -T 7 ^^2n 2 R 3 T r J du d 4 x ^9 9^' 9™' F% F; v , (3.6) 

Ub=0 

where we use the following definitions for the D7-brane tension T 7 and the trace over the repre- 
sentation matrices t a , 

T 7 =(2tt)VK)- 4 , (3-7) 
ti{t a t b ) =T r 5 ab . (3.8) 

In our case we have T r = 1/2. The overall factor 2% 2 R 3 comes from the integration over the 
5,6, 7-directions, which are the directions along the S 3 . 

Evaluating the DBI action given in (13.61) with the substitution rule (|3.5I) . we arrive at 

S D7 = -T 7 ^p^2n*R 3 T r 



x 



dud 4 * ^=99^9^ (4<V^ d W A h ~ 8^<Wr%^] A\, //) 



(3.9) 



where we use the short-hand notation [f = fx5 3c and neglect terms of higher than linear order in 
fx, and higher than quadratic order in A since both are small in our approach. 
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Up to the sum over flavor indices a, the first term in the bracket in (13.91 ) is reminiscent of the 
Abelian super-Maxwell action in five dimensions, considered already for the R charge current 
correlators in [2]. The new second term in our action arises from the non- Abelian nature of the 
gauge group, giving terms proportional to the gauge group's structure constants f abc in the field 
strength tensor F* v = 2d [fi A a u] + f abc A\ A%. 



D. Equations of motion 

We proceed by calculating the retarded Green functions for the action (13.91 ), following the 
prescription of Q3U as outlined in section 2 above. According to this prescription, as a first step we 
consider the equations of motion obtained from the action (|3.9I) . which are given by 



= 2d ll (y=gg^gr' <V^ 



+ / 



abc 



(3.10) 



It is useful to work in momentum space from now on. We therefore expand the bulk gauge fields 
in Fourier modes in the x l directions, 

A MA = f ^e--*° +ife -M>,£). (3.11) 

As we work in the gauge where A u = 0, we only have to take care of the components Ai with 

i = 0,1,2,3. 

For the sake of simplicity, we choose the momentum of the fluctuations to be along the x^ 
direction, so their momentum four- vector is k = (to, 0, 0, q). With this choice we have specified to 
gauge fields which only depend on the radial coordinate u, the time coordinate x and the spatial 
x 3 direction. 



1. Equations for Af - and A^-components 

Choosing the free Lorentz index in the equations of motion (13.101) to be v = a = 1, 2 gives two 
identical differential equations for A\ and A 2 , 

f tt,2 _ f Q 2 ft, ,.b 

= A a J' + Ljg + -P-Al + 2i—r bc 1 f-A c a , (3.12) 

/ U J U J 27Ti 

where we indicated the derivative with respect to u with a prime and have introduced the dimen- 
sionless quantities 

ro = , q = — — , m=-*— . (3.13) 

We now make use of the structure constants of SU(2), which are f abc = e abc , where e abc is the 
totally antisymmetric epsilon symbol with e 123 = 1. Writing out (13.121) for the three different 



9 



choices of a = 1, 2, 3 results in 



= Al" + LAI' + ^P-Al + 2^A\ , (3.15) 



= Al" + ^A 3 a ' + J q \ l . (3.16) 



a f a uf 

The first two of these equations are coupled, the third one is the same equation that was solved in 
the Abelian Super-Maxwell case 



2. Equations for Aq- and A^-components 

The remaining choices for the free Lorentz index v = 0, 3, u in (13.101 ) result in three equations 
which are not independent. The choices v = and v = u give 

n - A a " — — A a — —A a — n—f ahc -^—A c C\\l\ 
° " A ° uf A ° uf A3 \f* 2*T A *> (3 - 1?) 



= tt)< + qfAf + if abc ^A c > . (3.18) 



Solving (|3.18l) for Aq, differentiating it once with respect to u and using (13.171) results in equation 
(I3T01) for u = 3, 

= Af + + ^- 2 A% + + i\j ahc f= Al + 2i^- 2 r bc /- At . (3.19) 

j uj A uj z uj z znl uj A 2tt1 

We will make use of the equations (I3.17I) and (I3.18I) which look more concise. These equations of 
motion for Aq and A^ are coupled in Lorentz and flavor indices. To decouple them with respect 
to the Lorentz structure, we solve (13 .18b for A%' and insert the result into the differentiated version 
of (13.171) . This gives 

n _ A al " 4- A a " 4- ^ ~ ^ A af 4- 9? ro f abc ^ A cl n ?m 

0-A + uf A Q + uf2 A +2i uf2 f ^A,. (3.20) 

The equations for a = 1, 2 are still coupled with respect to their gauge structure. The case 
a = 3 was solved in [0]. We will solve (|3.20D for A^ and can obtain Af from (|3.18D . Note that it 
is sufficient for our purpose to obtain solutions for the derivatives of the fields. These contribute to 
(12.41) . while the functions A = f(u, k)A bdy (k) themselves simply contribute a factor of f(u, —k) 
which is one at the boundary. 



3. Solutions 

Generally, we follow the methods developed in [2], since our differential equations are very 
similar to the ones considered there. Additionally, we need to respect the flavor structure of the 
gauge fields. The equations for flavor index a = 3 resemble the ones analyzed in However, 
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those for a = 1,2 involve extra terms, which couple these equations. In our case the equations are 
coupled not only via their Lorentz indices, but also with respect to the flavor indices. We already 
decoupled the Lorentz structure in the previous section. As shown below, the equations of motion 
which involve different gauge components will decouple if we consider the variables 

Xi= Al + iAl 

(3.21) 

X = Aj-iAl 

Here the A], A\ are the generally complex gauge field components in momentum space. Note 
that up to SU(2) transformations, the combinations (13.211) are the only ones which decouple the 
equations of motion for a = 1,2. These combinations are reminiscent of the non-Abelian SU(2) 
gauge field in position space, 

A-A^-U A? A "~ iA A (3 22) 

A ~ Ai 2 "2 \Al + %A\ -A\ ) ■ (3 - 22) 

The equations of motion for the flavor index a = 3 were solved in [0]. To solve the equations 
of motion for the fields A" with a = 1, 2, we rewrite them in terms of X { and X { . Applying the 
transformation (13.211) to the equations of motion (13.141) and (13.151) and the a = 1, 2 versions of 
(ET201) and (I3T81) leads to 



f x >, ft) 2 ~ /q 2 = F 2mro 
f a uP 
juf)' Y „ tu 2 - /q 2 = F 2mtu 
uf 0+ uP 
= (tv T m)X^ + qfX> 3 , (3.25) 



K + hK + ~ J "J — X a , a = 1, 2, (3.23) 

= X"' + K -^-X'> + ' U ~ / V" tmu ^, (3.24) 



where the upper signs correspond to X and the lower ones to X. 

We see that some coefficients of these functions are divergent at the horizon u = 1. Such 
differential equations with singular coefficients are generically solved by an ansatz 

X i = (l-ufF(u), X i = (l-ufF(u), (3.26) 

with regular functions F(u) and F(u). To cancel the singular behaviour of the coefficients, we 
have to find the adequate (3 and (3, the so-called indices, given by equations known as the indicial 
equations for (3 and (3. We eventually get for all Xj and X{ 



, itt) L 2m ~ , iro I 2m 
P = ±—Jl-— p = ± Jl+ (3.27) 

2 V tf 2 V tf 

Note that these exponents differ from those of the Abelian Super-maxwell theory [2] by a de- 
pendence on v^fo in the limit of small frequencies (ro < m). In the limit of vanishing chemical po- 
tential m — > 0, the indices given in [|20 are reproduced from (13.271 ). In order to solve (13.231 ), (13.241 ) 
and (13.251) . we wish to introduce a series expansion ansatz in the momentum variables ro and q. In 
fact, the physical motivation behind this expansion is that we aim for thermodynamical quantities 
which are known from statistical mechanics in the hydrodynamic limit of small four-momentum 
k. So the standard choice would be 



F(u) =F + xvF 1 + q 2 d + 



(3.28) 
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On the other hand, we realize that our indices will appear linearly (and quadratically) in the differ- 
ential equations' coefficients after inserting (13.261) into (13.231) . (13.241) and (|3.25l) . The square root 
in (3 and (3 mixes different orders of ro. In order to sort coefficients in our series ansatz, we assume 
to < m and keep only the leading to contributions to (3 and (3, such that 



/torn ~ , . /torn 

, f}~^—. (3.29) 

This introduces an additional order O(vo l l 2 ), which we include in our ansatz (13.281) yielding 

F(u) =F + to 1/2 F 1/2 + toF x + q 2 G x + . . . , (3.30) 

and analogously for the tilded quantities. If we had not included O(vo l l 2 ) the resulting system 
would be overdetermined. The results we obtain by using the approximations (13.291) and (|3.30l) 
have been checked against the numerical solution for exact (3 with exact F(u). These approxima- 
tions are useful for fluctuations with q, to < 1 (see subsection IB 41 in Appendix |B]). Note that by 
dropping the 1 in (13.271) we also drop the Abelian limit. 

Consider the indices (13.271 ) for positive frequency first. In order to meet the incoming wave 
boundary condition introduced in section [TTl we restrict the solution (3 to the negative sign only. 
For the approximate (3 in (13.291) we therefore choose the lower (negative) sign. This exponent 
describes a mode that travels into the horizon of the black hole. In case of (3 we demand the mode 
to decay towards the horizon, choosing the lower (positive) sign in (13.291 ) consistently. Note that 
for negative frequencies uj < the indices (3 and (3 exchange their roles. 

Using (13.291 ) in (13.261) and inserting the ansatz into the equations of motion, we find equations 
for each order in q 2 and vo separately. After solving the equations of motion for the coefficient 
functions F , F 1 / 2 , F\ and G\, we eventually can assemble the solutions to the equations of motion 
for X as defined in (13.211) . 

X(u) = (1 - uf F{u) = (1 - uf (F Q + VnF 1/2 + nF 1 + q 2 d + ...). (3.31) 

and a corresponding formula for X(u) from the ansatz (13.261) . 

Illustrating the method, we now write down the equations of motion order by order for the 
function X a . To do so, we use (13.311) with (13.291) in (13.231) with the upper sign for X a . Then we 
examine the result order by order in vo and q 2 , 

O(const): = ^' + ^, (3.32) 



O(VS) : = F{> /2 + f jF[ /2 - ^F,; - y|iF , (3.33) 



OM: 0^ + ^-^-^-«iz^ Fo , (3.34) 



0(q 2 ): Q = G'l+ f ^G' 1 -—F. (3.35) 

At this point we observe that the differential equations we have to solve for each order are 
shifted with respect to the solutions found in [0]. The contributions of order ro n in [0] now show 
up in order ro"/ 2 . Their solutions will exhibit factors of order /W 2 . 
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Solving the system (13.321) to (13.351) of coupled differential equations is straightforward in the 
way that they can be reduced to several uncoupled first order ordinary differential equations in 
the following way. Note that there obviously is a constant solution F = C for the first equation. 
Inserting it into (13.331 ) and (13.351 ) leaves us with ordinary differential equations for F[ , 2 and G[ 
respectively. Using the solutions of F and Fi/ 2 in (|3.34l) gives one more such equation for F[. 

To fix the boundary values of the solutions just mentioned, we demand the value of F{uh = 1) 
to be given by the constant F and therefore choose the other component functions' solutions 
such that lim n ^i Fi/ 2 = 0, and the same for F\ and G\, The remaining integration constant C is 
determined by taking the boundary limit u — ► of the explicit solution (|3.31l) . making use of the 
second boundary condition 

liml(u) = X bdy , (3.36) 

see appendix |B] Eventually, we end up with all the ingredients needed to construct the gauge 
field's fluctuations X(u) as in (|3.31l) . 

We solve the equations (13.231) with lower sign for X a and (|3.24l) for X' and its tilded partner in 
exactly the same way as just outlined, only some coefficients of these differential equations differ. 
The solution for X' 3 is then obtained from (13.251) . 

All solutions are given explicitly in Appendix |B] together with all other information needed to 
construct the functions X a , X a , X' Q , X' Q , X' 3 and X 3 . 



IV. ISOSPIN DIFFUSION AND CORRELATION FUNCTIONS 



A. Current correlators 



In this section we obtain the momentum space correlation functions for the gauge field compo- 
nent combinations X and X defined in equation (13.211) . Recall that the imaginary part of the 
retarded correlators essentially gives the thermal spectral functions (see also section [II]). The 
following discussion of the correlators' properties is therefore equivalent to a discussion of the 
corresponding spectral functions. 

First note that the on-shell action gets new contributions from the non-Abelian structure, 



D7 



(2W) 2 2 o 
T 7 K . 27i 2 R 3 T r 



x 2 



2 



(4.1) 



u b =0 




-99°°9 33 A^A b 0] 



where j, j' = 0,1,2,3 and the index u denotes the radial AdS -direction. Up to the sum over 
flavor indices, the first term in the bracket is similar to the Abelian Super-Maxwell action of 
The second term is a new contribution depending on the isospin chemical potential. It is a contact 
term which we will neglect. The correlation functions however get a structure that is different 
from the Abelian case. This is due to the appearance of the chemical potential in the equations of 
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motion and their solutions. Writing (14.11) as a function of X and X results in 



Si 



D7 



-T 7 ^2n^T r 



x 2 



2 
d 4 g 

(2^ 



-(X/X^ + XjX^+^f, 



M b =0 



(4.2) 



4g/i J du^g w g M ( X [0 X 3] + X [3 X 0] 
o 



In order to find the current correlators, we apply the method outlined in section [XT] to (14.21) . 
with the solutions for the fields given in appendix [Bj As an example, we derive the correlators 
Gqo — (Jo(q)Jo(—<T) ) an d Gq = (Jo(q) Jo(—q)) of the flavor current time components J and J , 
coupling to the bulk fields X and X , respectively. Correlation functions of all other components 
are derived analogously. For the notation see appendix lAl 



1. Green functions: Calculation 
First, we extract the factor B(u) of (12.21) . 

B(u) = 



(27ra') 2 2 3 , — - 00 
i 7 27r ^ l r y/-gg g . 



(4.3) 



The second step, finding the solutions to the mode equations of motion, has already been per- 
formed in section UlID 31 In the example at hand we need the solutions X and X . From (13.311 ) 
and from appendix |B] we obtain 



X ' 



u, 



q 2 X My + tt,qX 3 b 



V2rnfo + rom In 2 + q 2 



m 



K>— I tx z + 31n z 2 + 31n z (l + u) + 6 In 2 In 



u 



X ' 



+12Li 2 (l -u) + 12Li 2 (-u 
fl-uV 



12Li 2 



q 2 X bdy + raqX 3 bdy 



1 — u 



1 + u 
+ q 2 ln 



2v^_ 
u+ 1 



u + 1 



(4.4) 



2u 



v / 2mfo + rom In 2 - q 2 



2u 2 



u 
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Ti 



+ 31n 2 2 + 31n 2 (l + w) + 61n21n 



u 



-12Li 2 (l -u) + 12Li 2 (-u) 



12Li 2 



1 — u 



1 + u 

2 



(4.5) 



q^ln 



u + 1 
2u 



Note that we need the derivatives to apply (I2.4I ). 

Now we perform the third step and insert (14.31) , (14.41 ) and (14.51 ) into (I2.4I ). Our solutions X 
and X replace the solution f(u, k) and f(u, —k) in (12.41) . The resulting expression is evaluated 
at Ub = 0, which comes from the lower limit of the w-integral in the on-shell action (14.21) . At small 
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it = e — > 0, (l4~4l) and (1431) give 

lim x ' = - q^o dy + _ / 2 ~ bdy ~ bdy \ 

lim X ' = q^g dy + ^qXg _ + Hm / + ^ g ^ 

zv2mtt) + ram In 2 - q 2 V / 

In the next to leading order of (14.61) and (14.71) there appear singularities, just like in the Abelian 
Super-Maxwell calculation [0, equation (5.15)]. However, in the hydrodynamic limit, we consider 
only the finite leading order. 



2. Green functions: Results 

Putting everything together, for the two Green functions for the field components X , X given 
indUDby 

X = Al + iA 2 , X = Al-iAl 

we obtain 

™f .. . (4.8) 



57T 



V2mra — q 2 — ram ln2 



N C T 2nT q 2 

^ iv / 2mra - q 2 + ram ln2 



Goo = ^ . ;== V. • (4-9) 



These are the Green functions for the time components in Minkowski space, perpendicular to the 
chemical potential in flavor space. All Green functions are obtained considering hydrodynamic 
approximations in C(ra 1/ ' 2 , ra, q 2 ), neglecting mixed and higher orders 0(ra 3//2 , ra 1/,2 q 2 , q 4 ). 

The prefactor in (14.81) . (14.91) is obtained using T7 as in (13.71) . and carefully inserting all metric 
factors, together with the standard AdS/CFT relation R A = A-Kg s N c a' 2 . As in other settings with 
flavor Q26Q, we concordantly get an overall factor of N c , and not iV 2 , for all correlators. Contrary 
to those approaches, we do not get a factor of Nj when summing over the different flavors. This is 
due to the fact that in our setup, the individual flavors yield distinct contributions. Most striking is 
the non-trivial dependence on the (dimensionless) chemical potential m in both correlators. Note 
also the distinct structures in the denominators. The first one (14.81 ) has an explicit relative factor 
of i between the terms in the denominator. In the second correlator (14.91) there is no explicit factor 
of i. The correlator (14.81 ) has a complex pole structure for to > 0, but is entirely real for to < 0. 
On the other hand, (14.91 ) is real for to > but develops a diffusion structure for uj < 0. So the 
correlators G q and Gq essentially exchange their roles as uj changes sign (see also Fig. Q]). We 
find a similar behavior for all correlators G^and G~- { with j, I — 0, 1, 2, 3. This behavior is a 

consequence of the insertion of (^(ra 1 / 2 ) in the hydrodynamic expansion (13.301) . 

We assume m to be small enough in order to neglect the denominator term of order C(ram) 
(9(v / ram, q 2 ). Moreover, using the definitions of ra, q and m from (13.131) we may write (|4.8|) 
and (1491 as 



G ~ = - 7 ^ 7 = ; (4.10) 



N C T q 2 ^J 
8iTy/2jl uj + q 2 D(uj) 
N C T q 2 ^J 
8n^/2ji iuj — q 2 D(uj) 
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where the frequency-dependent diffusion coefficient D{oj) is given by 

D < w) = ^Sr- (4 - 12) 

We observe that this coefficient also depends on the inverse square root of the chemical potential /x. 
Its physical interpretation is discussed below in section HVB I 
In the same way we derive the other correlation functions 



Sttv^I Q(fj,q) 33 Stt^ Q(u,q) 



G = ^ y/uq(u-n) _ N C T ^q(uj + fi) 



iV c T u^q _ N C T u^q 

30 %-K^2jI Q^^q)' 30 STTv^g^g)' 

with the short-hand notation 

Q(u,q) = ioo - q 2 D(u) } Q(u, q) = oo + q 2 D(u) . (4.16) 

Note that most of these functions are proportional to powers of q and therefore vanish in the limit 
of vanishing spatial momentum q — > 0. Only the 33-combinations from (14.131) survive this limit. 
In contrast to the Abelian Super-Maxwell correlators [2] given in appendixO it stands out that our 
results (14.101 ), (14.1 II) and (14.131 ) and (14.151 ) have a new zero at cu = /x or — /x. Nevertheless bear in 
mind that we took the limit cu < \i in order to obtain our solutions. Therefore the apparent zeros at 
±/x lie outside of the range considered. Compared to the Abelian case there is an additional factor 
of \fu. The dependence on temperature remains linear. 

In the remaining X-correlators we do not find any pole structure to order y/u, subtracting an 
0(q 2 ) contribution as in [|20, 

G ri =G^ = ^^^, (4.17) 

G~n =G~ 22 = --^^JjEj. (4.18) 

As seen from (14.171) . G a & (with a = 1, 2) are real for negative cu and imaginary for positive to. The 
opposite is true for G^a, as is obvious from the relative factor of i. 

The correlators of components, pointing along the isospin potential in flavor space (a = 3), are 
found to be 

Ga « a ° = iu> - zw ' Ga1a% = Ga1a1 = ^T™ - zw ' (4 ' 19) 

G AlAl = G A%A% = , G AUl = -j- - , (4.20) 

with the diffusion constant D = l/(2nT) . Note that these latter correlators have the same 
structure but differ by a factor 4/iV c from those found in the Abelian super-Maxwell case [0] (see 
also (IC1I) and (IC2I) ). In particular the correlators in equation (14.191 ) do not depend on the chemical 
potential. 
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To analyze the novel structures appearing in the other correlators, we explore their real and 
imaginary parts as well as the interrelations among them, 

NT a 2 

ReG o5 (u; > 0) = ReG 5o (cj < 0) = - ^ , (4.21) 



NT o 4 

ReG o5 (^<0)= ReG~ 00 (u>0) = - J,,.,.,,., , (4.22) 



8vr ^/2^\J\+q^/(2-KT)' 
N C T q A 
16vr 2 2\x + g 4 /(2vrT) 2 



iV T 

lmG ~ (u<0)= -lmG~ 00 (u>0)= ^ __^£1^L_ , (4 .23) 
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V2A*I 


a; 




2/i 


a; 


+ g 4 /(27rT) 2 



and 



ImG o5 (w>0) = 0, ImG 5o (w<0) = 0. (4.24) 

Now we see why, as discussed below (14.1 II) . G q and Gq exchange their roles when crossing the 
origin at uj = 0. This is due to the fact that the real parts of all G.j-and G~- x are mirror images 
of each other by reflection about the vertical axis at to = 0. In contrast, the imaginary parts are 
inverted into each other at the origin. Figured] shows the real and imaginary parts of correlators 
G q and Gg . The different curves correspond to distinct values of the chemical potential /i. The 
real part shows a deformed resonance behavior. The imaginary part has a deformed interference 
shape with vanishing value for negative frequencies. All curves are continuous and finite at uj = 0. 
However due to the square root dependence, they are not differentiable at the origin. Parts of the 
correlator which are real for positive uj are shifted into the imaginary part by the change of sign 
when crossing to = 0, and vice versa. 

To obtain physically meaningful correlators, we follow a procedure which generalizes the 
Abelian approach of [6]. In the Abelian case, gauge-invariant components of the field strength ten- 
sor, such as E a = ojA a , are considered as physical variables. This procedure cannot be transferred 
directly to the non-Abelian case. Instead, we consider the non-local part of the gauge invariant 
trF 2 which contributes to the on- shell action (14.11) . In this action, the contribution involving the 
non-Abelian structure constant - as well as /i - is a local contact term. The non-local contribution 
however generates the Green function combination 

Ga 1 a 1 + G A 2 A 2 + G A 3 A 3 . (4.25) 

We take this sum as our physical Green function. This choice is supported further by the fact 
that it may be written in terms of the linear combinations (13.211 ) which decouple the equations of 
motion. For example, for the time component, written in the variables X , X given by (13.211) . 
the combination (14.251) reads (compare to (14.21) ) 

G od + G qq + G A 3 A 3 . (4.26) 



The contribution from £^3^3 is of order 0(/jP), while the combination for the first two flavor 
directions, G q + Gg , is of order 

We proceed by discussing the physical behavior of the Green function combinations introduced 
above. G A 3 A s is plotted in Fig. [2] on the right. Its frequency dependence is of the same form as in 
the Abelian correlator obtained in [2], as can be seen from (IC1I) . Since we are interested in effects 
of order O(fx), we drop the third flavor direction a = 3 from the sum (|4.26l) in the following. It 
is reassuring to observe that the flavor directions a = 1, 2 which are orthogonal to the chemical 
potential, combine to give a correlator spectrum qualitatively similar to the one found in [2] for 
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FIG. 1: Real (left plot) and imaginary part (right plot) of the correlator G Q q as a function of frequency uj/T 
at different chemical potential values fi/T = 0.5 (solid line), /i/T = 0.3 (short-dashed line) and fx/T = 
0.1 (long-dashed line). The corresponding plots for the correlator G Q q would look like the mirror image of 
the ones given. The real part would be reflected about the vertical axis at u> = 0, the imaginary part would 
be reflected about the origin. All dimensionful quantities are given in units of temperature. The numerical 
values used for the parameters are q/T = 0.1, N c = 100. 

the Abelian Super-Maxwell action (see Fig. [2]). However, we discover intriguing new effects such 
as the highly increased steepness of the curves near the origin due to the square root dependence 
and a kink at the origin. 

We observe a narrowing of the inverse resonance peak compared to the form found for the 
Abelian Super-Maxwell action (and also compared to the form of our G A s A s, as is seen from 
comparing the left with the right plot in Fig. [2]). At the origin the real and imaginary part are 
finite and continuous, but they are not continuously differentiable. However, the imaginary part 
of G A 3 A 3 has finite steepness at the origin. The real part though has vanishing derivative at uj = 0. 
Note that the imaginary part of flavor directions a = 1, 2 on the left plot in Fig. |2] never drops 
below the real part. In the third flavor direction, as well as in the Abelian solution, such a drop 
occurs on the positive oaxis. 

The correlators G 33 , G 33 , G 03 ~ and Gq 3 have the same interrelations between their respective 
real and imaginary parts as G q and Gq . Nevertheless, their dependence on the frequency and 
momentum is different, as can be seen from (14.131) to (14.151) . A list of the 33-direction Green 
functions split into real and imaginary parts can be found in appendix iDl 

Thermal spectral functions in different directions are compared graphically in appendix [0 

B. Isospin diffusion 

The attenuated poles in hydrodynamic correlation functions have specific meanings (for ex- 
emplary discussions of this in the AdS/CFT setup see e.g. [7J], jSOQ ). In our case we observe an 
attenuated pole in the sum G ^ + Gq at uj = 0. As can be seen from the plots in Fig. |2l our pole 
lies at Recu = 0. This structure appears in hydrodynamics as the signature of a diffusion pole 
located at purely imaginary uj. Its location on the imaginary cu-axis is given by the zeros of the 
denominators of our correlators as (neglecting 0(u> , q A )) 



(4.27) 
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FIG. 2: In the left plot the sum of both correlators in 00-directions is split into its imaginary (dashed 
line) and real (solid line) part and plotted against frequency. For comparison the right plot shows the 
corresponding real and imaginary parts for the G A 3 A 3. It is qualitatively similar to the Abelian correlator in 
i = Lorentz direction computed from the Super-Maxwell action in |2fl. Note the different frequency scales 
in the two plots. The curves in a = 1, 2-directions are much narrower due to their square root dependence 
on lo. Furthermore they have a much larger maximum amplitude. All dimensionful quantities are given in 
units of temperature. The numerical values used for the parameters are, as in Fig.Q] q/T = 0.1, N c = 100 
and only in the left plot fi/T = 0.2. 



Squaring both sides of (|4.27l) we see that this effect is of order 0(q A ). On the other hand, looking 
for poles in the correlator involving the third flavor direction G A z A z, we obtain dominant contri- 
butions of order 0(q 2 ) and 0(/jP) (neglecting 0(uj 2 , q 2 )) 

»--*sj^ («*> 

This diffusion pole is reminiscent of the result of the Abelian result of |Q] given in appendix dCT). 
As discussed in section HV A[ we consider gauge-invariant combinations G q + Gq + G A 3 A 3. In 
order to inspect the non- Abelian effects of order O(fi) showing up in the first two correlators in 
this sum, we again drop the third flavor direction which is of order O(fjP). 

Motivated by the diffusion pole behavior of our correlators in flavor-directions a = 1, 2 corre- 
sponding to the combinations X, X (see (14.271) ), we wish to regain the structure of the diffusion 
equation given in (12.91) , which in our coordinates (k = (u> , , , q )) reads 

iuj J = D(uj)q 2 J Q . {A 29) 

Our goal is to rewrite (14.271 ) such that a term 0{uS) and one term in order 0(q 2 ) appears. Further- 
more there should be a relative factor of —i between these two terms. The obvious manipulation 
to meet these requirements is to multiply (14.271) by y/uj in order to get 



^-^2^m- (430) 



Comparing the gravity result (14.301) with the hydrodynamic equation (14.291) . we obtain the 
frequency-dependent diffusion coefficient 

D <"> = iSsrT- <4 ' 31) 
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Our argument is thus summarized as follows: Given the isospin chemical potential as in (11.11) . 
(13.51) . Jo from ( 12.91) is the isospin charge density in (14.291) . According to (14.291 ). the coeffi- 
cient (14.311) describes the diffusive response of the quark-gluon plasma to a gradient in the isospin 
charge distribution. For this reason we interprete D{u) as the isospin diffusion coefficient. 

Near the pole, the strongly-coupled plasma behaves analogously to a diffractive medium with 
anomalous dispersion in optics. In the presence of the isospin chemical potential, the propaga- 
tion of non-Abelian gauge fields in the black hole background depends on the square root of the 
frequency. In the dual gauge theory, this corresponds to a non-exponential decay of isospin fluc- 
tuations with time. 

The square root dependence of our diffusion coefficient is valid for small frequencies. As 
long as u/T < 1/4, the square root is larger than its argument and at u/T = 1/4, the differ- 
ence to a linear dependence on frequency is maximal. Therefore in the regime of small frequen- 
cies u/T < 1/4, which is accessible to our approximation, diffusion of modes close to 1/4 is 
enhanced compared to modes with frequencies close to zero. 



V. CONCLUSION 



In this paper we have considered a relatively simple gauge/gravity dual model for a finite tem- 
perature field theory, consisting of an isospin chemical potential fj, obtained from a time component 
vev for the SU(2) gauge field on two coincident brane probes. We have considered the D7-brane 
embedding corresponding to vanishing quark mass, for which fi is a constant, independent in 
particular of the radial holographic coordinate. The main result of this paper is that this model, 
despite its simplicity, leads to a hydrodynamical behavior of the dual field theory which goes be- 
yond linear response theory. We find in particular a frequency-dependent diffusion coefficient 
with a non-analytical behavior. Frequency-dependent diffusion is a well-known phenomenon in 
condensed matter physics. Here it originates simply from the fact that due to the non-Abelian 
structure of the gauge field on the brane probe, the chemical potential replaces a time derivative in 
the action and in the equations of motion from which the Green functions are obtained. 

Of course the calculation presented has some limitations as far as the approximations made are 
concerned. This applies in particular to the approximation (13.291 ) of the so-called indices in the 
ansatz for solving the equations of motion. Here we have dropped the constant present under the 
square root and used an expression proportional to the square root of the frequency. This allows 
for a closed solution without having to use numerics. However using this approximation we have 
dropped the Abelian limit. This leads ultimately to the square root dependence of the diffusion 
coefficient on the frequency. This dependence is unphysical for u — ► 0, since the diffusion coeffi- 
cient is expected to be non-zero for zero frequency. We expect physical behavior to be restored if 
the Abelian limit is included in the calculation. To avoid the approximation described, this requires 
a numerical approach. A suitable solution method for all momenta and all frequencies has been 
presented in [7] and in Il49ll . We are going to study the application of this method to the model 
presented here in the future. 
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APPENDIX A: NOTATION 



The five-dimensional AdS Schwarzschild black hole space in which we work is endowed with 
a metric of signature (—,+,+,+,+), as given explicitly in (13.31) . We make use of the Einstein 
notation to indicate sums over Lorentz indices, and additionally simply sum over non-Lorentz 
indices, such as gauge group indices, whenever they occur twice in a term. 

To distinguish between vectors in different dimensions of the AdS space, we use bold symbols 
like q for vectors in the the three spatial dimensions which do not live along the radial AdS 
coordinate. Four-vectors which do not have components along the radial AdS coordinates are 
denoted by symbols with an arrow on top, as q. 

The Green functions G = (J J) considered give correlations between currents J and J. These 
currents couple to fields A and A respectively. In our notation we use symbols such as G A a A b to 
denote correlators of currents coupling to fields A\ and A\, with flavor indices a, b and Lorentz 
indices k, I = 0, 1, 2, 3. For the gauge field combinations X k and Xi given in (13.211) , we obtain 
Green functions G k j denoting correlators of the corresponding currents. 

APPENDIX B: SOLUTIONS TO EQUATIONS OF MOTION 

Here we explicitly write down the component functions used to construct the solutions to the 
equations of motion for the gauge field fluctuations up to order ro and q 2 . The functions themselves 
are then composed as in (13.311) . 

The solutions for the components with flavor index a = 3 where obtained in 



1. Solutions for X a , X a and A^ a 



The function X a {u) solves (13.231) with the upper sign and is constructed as in (13.311 ) from the 
following component functions, 
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where the constant C can be expressed it in terms of the field's boundary value X bdy 
\im u ^ Q X(u, k), 



u, K i /mtu fix 2 In 2 

C = X bdy [ 1 + J — l n 2 + mw (y + — 



)+^q 2 + 0(tt, 3 / 2 ,q 4 )j . (B6) 



The solutions of the equations of motion (13.231) with lower sign for the functions X a (u) are 
given by 
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with C given by 



C = X hdy 1 - 



ln2-mm(^ + ^+^q 2 + 0(m 3 / 2 ,q 4 )') , (B12) 



so that lim u ^o X(u, k) = X bdy . 

The solution for A^ solves (13.161) up to order m and q 2 with boundary value (A 3 ) y . It is 
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2. Solutions for X' and Ag' 



Here we state the solutions to (13.201) . This formula describes three equations, differing in 
the choice of a = 1, 2, 3. The cases a = 1, 2 give coupled equations which are decoupled by 

1 2 " — ' 

transformation from A ' to X and X . The choice a = 3 gives a single equation. 
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The function X' Q is solution to (13.241) with upper sign. We specify the component functions as 
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where the constant C can be expressed in terms of the field's boundary value X bdy 
lim n ^ ^(«, k), 
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V2mro + mro In 2 + q 2 
To get the function X' Q , we solve (13.241) with the lower sign and obtain 







to m ^ , . 
->\l — + 0{w) 



F = C, 

~ ~ fm 2u 2 

Fi/o = iC \ — In , 

' V 2 1 + u' 



F 1 = C™ 
12 



7r 2 + 31n 2 2 + 31n 2 (l + u) + 61n21n 



d = Cln 



+ 12 I Li 2 (l -u) + Li 2 (-u) - Li 2 
1 + u 



1 + u 

1 — It 



2m 



(B20) 

(B21) 

(B22) 
(B23) 

(B24) 

(B25) 
(B26) 



where the constant C can be expressed it in terms of the field's boundary value X My 

\im u ^ X(u, k), 



C 



q 2 X bdy + raqX 3 b 



iv / 2mro + mtt) In 2 — q 2 
The solution for (13.201) with a = 3 is the function Aq, given by 



Af = (l-u)~ 



qM b dy + tvqAf* 



itv — c\ z 



in> . 2-u 2 



1 + — In ■ 



2 1 + u 



+ q 2 ln 



1 + u 
2u 
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FIG. 3: These plots show the real and imaginary part of the function F(u) which is part of X a = (1 — 
UfF(u). The solid line depicts the analytical approximation, obtained in this paper. As a check we solved 
the equations of motion for F(u) numerically. They are drawn as dashed lines. In this example we used 
T = 1. The numerical solution was chosen to agree with the analytical one at the horizon and boundary. 

3. Solutions for XL X'n and A 3 ,' 



We give the derivatives of X 3 and X 3 as 



X' 



ro — m 
ro + m 



The solution for A|' is 



At 



ro 



A 3 ' 



o • 



(B29) 
(B30) 

(B31) 



4. Comparison of numerical and analytical results 

As an example, in Fig. [3] we show the numerical and analytical solutions for the function F(u) 
in X a = (1 — uYF{u). Here we compare the numerical result for F{u) obtained from the ansatz 
(13.261) with (13.271) in (13.231) with the analytically obtained approximation given above in (IB 1 1) to 
(IB61) . 



APPENDIX C: ABELIAN CORRELATORS 

For convenient reference we quote here the correlation functions of the Abelian super-Maxwell 
theory found in [|2|]. The authors start from a 5-dimensional supergravity action and not from a 
Dirac-Born-Infeld action as we do. Therefore there is generally a difference by a factor N c /A. 
Note also that here all Nf flavors contribute equally. In our notation 

iN 2 Ti,> fr ab N 2 Tn 2 h ab 

rab _ rab N 2 Tuq 5 ab b N 2 Tu 2 5' 

G ° 3 " Gso " 16tt(^ - Dq 2 ) + " ' ' G33 " 16n(tu - Dq" + " ' ' ' ' 



where D = 1/(2ttT) 
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APPENDIX D: CORRELATION FUNCTIONS 



In this appendix we list the real and imaginary parts of the flavor currents in the first two 
flavor-directions a = 1, 2 and in the third Lorentz-direction coupling to the supergravity-fields X$ 
and X 3 (as defined in (13.211) ). 



Re{G 3 3^ > 0)} 
Im{G 3 3^ > 0)} 
Re{G 33 (^ < 0)} 



and 



Re{G 3 > < 0)} 
-Im{G 3 > < 0)} 
Re{G~ 33 (u > 0)} 



N c q 2 (u 2 + fi\uj\) 
16tt 2 [2fi\u\ +q 4 /(2irT) 2 } 



N C T 


LO 


(oo 2 + n 


lv\) 


8vr [2(1 


LO 


+ g 4 /(2vrT) 2 ] 



N C T( 


UJ 2 -(1 




8tt 


y/2fl 


LO 


+ q 2 /(2nT) 



lm{G^(co < 0)} = 0, Im{G3 3 (cc; > 0)} = 0. 



(Dl) 
(D2) 
(D3) 

(D4) 



APPENDIX E: THERMAL SPECTRAL FUNCTIONS 



We include here a comparision of the sizes of spectral functions in distinct flavor- and Lorentz- 
directions (see also (12.11) in section [III). 




FIG. 4: Here the thermal spectral functions in distinct Lorentz- and flavor-directions are plotted against 
frequency to/T in units of temperature. In the left plot the chemical potential was chosen to be (i/T = 0.7, 
in the right one fi/T = 0.2. Flavor-directions a = 1, 2 are summed and displayed as one curve. The 
frequency-dependence of 00- (solid red line) and 03-Lorentz-directions (short-dashed blue line) is shown. 
By the dotted line we denote the spectral curve in 11- or 22-directions. This curve was scaled by a factor 100 
in order to make it visible in these plots. The third flavor-direction is only plotted for the spectral function in 
Lorentz-directions 00 (long-dashed green curve). We do not show the 33-direction spectral function which 
has a square root dependence and is comparable in size with the 11-direction. 
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